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Powerful extended radio galaxies like Cygnus A can be used for cosmological tests. 
The characteristic, time-independent length _D* for several radio sources is calculated 
and compared to the average physical length for a sample of radio galaxies. The 
ratio of these two lengths as a function of redshift is used to constrain cosmological 
parameters. Cygnus A is the only low- redshift (z < 0.5) radio galaxy for which we 
have an estimate for D*. We comment on the sensitivity to this one low-redshift point, 
and results including and excluding Cygnus A are shown. A low density universe is 
favored, and the one free model parameter is relatively well constrained. 



1. Introduction 

The use of the physical sizes of radio sources as a cosmological test has been 
considered by many [1,2,3,4] since radio sources are observed to fairly high red- 
shift. The coordinate distance (a r), used to calculate the size of a source is 
a function of the density parameter = 8irGp /3H , and the normalized 
cosmological constant A = A/3H , where p is the mean mass density of the 
universe, H is the Hubble constant, and A is the cosmological constant. It has 
been shown that radio sources are not, strictly speaking, standard yardsticks [4]. 
However, if the evolution with redshift of physical sizes could be predicted and 
compared to observations, then radio sources could be used for cosmography [5]. 

Daly [6] presents a method where a characteristic size D* (called /* by Daly 
[6,7]) for a powerful extended radio source can be calculated from observables 
assuming a value for one model parameter f3. The evolution of these chracteristic 
sizes for a chosen sample should, for the correct parameters , A ), track 
the evolution of the average physical size of a large sample of similar sources 
(see §2). The cosmological dependence of the ratio of these sizes goes as a power 
greater than one of (a r) for realistic values of f3 (see §4). Thus, as we compare 
samples at high redshift, the dependence of (a r) on and A will allow us to 
constrain the cosmological parameters. It should be noted that the results of 
this analysis are insensitive to the value of H (see §4); here H is parameterized 
in the usual way: H = 100 h km s" 1 Mpc" 1 . 

To calculate D* for a radio source, one must have the lobe propogation velocity 
vl } the lobe magnetic field strength B } and the lobe width (see §2). These 
parameters can be extracted from multi-frequency observations of radio sources 
[8,9]. 



2 E. J. Guerra & R. A. Daly: Cygnus A and Other SCR Sources 



2. The Characteristic Size D* and the True Physical Size D 

Daly [6] describes how to calculate a characteristic, time-independent length 
scale for a radio source. If a radio source has constant lobe propogation ve- 
locity vl with a stable direction of outflow, and is active for a lifetime £*, the 
characteristic length is 

D*=v L t*. (2.1) 

One can think of D* as the core-lobe size at the end of the source's lifetime, or as 
the average lobe-lobe separation. This method may be applied to sources that 
are supersonic propagators with no significant backflow or bridge distortions, 
such as the high power FRII radio galaxies (see §3 and [6]). Note that present 
observations are consistent with vl being roughly constant over the lifetime of 
a particular source [7,9,10]. 

If the luminosity in directed kinetic energy, Lj, is roughly constant over the 
lifetime of the source, then ~ Ei/Lj where E{ is the initial energy of the 
central engine that powers the jets. Note that E{ and Lj will depend only on 
processes occuring near or at the central engine. A power law is assumed between 

and Lj 

U oc L~^\ (2.2) 

which is equivalent to assuming E{ oc Lj 1- ' 9 / 3 (see [6]). The determination of f3 
has implications for models of energy extraction from the central engine [6,7,11]. 
Combining eqs. (2.1) and (2.2), it is clear that 

D* oc vlLj~^ 3 . (2.3) 

The radio lobe propagates supersonically and does P dV work on the ambient 
medium [12]. The pressure P of the lobe is given by the strong shock jump 
conditions P = \p a "VL 2 where p a is the mass density of the ambient medium. 
The energy input to the hotspot and lobe in a time interval St is SE ~ Lj St. 
Equating SE to P dV with volume element dV = rra^viSt gives a relation 
between vl and Lj: 

Lj oc n a a L 2 v L 3 oc (2.4) 



where n a is the number density of the ambient media and 
Combining eqs. (2.3) and (2.4) 



k = (n a a L 2 ) 1/3 . (2.5) 



£>* oc ifc'V 1 "' 9 . (2.6) 

For supersonic flow, the lobe is ram pressure confined [13], so the pressure of 
the lobe P oc n a VL 2 . The lobe pressure will be oc B 2 if the magnetic field B and 
the relativistic electrons are near equipartition. Thus, the ambient gas density 
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can be estimated as [6] 



n a oc (B/v L ) 



(2.7) 



Combining eqs. (2. 5), (2. 6), and (2.7), gives an expression for D* in terms of 
observables. 



For our purposes, the normalization of D* is unimportant and arbitrary. Note 
that in theory, D* is time-independent for a given source and thus is the same 
no matter when the source is observed during its lifetime. 

The core-lobe size of a radio source D is given by D = v^t where t is the age 
of the source and vl is assumed to be constant. Thus for a given source D = 
(t/t*)D*. For the sources we consider, the lobe-lobe size is well approximated 
by 2D (see §4). 

In a given redshift range, one expects the D* } s for a sample to be peaked 
around some central value. Taking the D* } s in this redshift range to be approx- 
imately equal to a constant, D Z} it follows that 



If we assume the distribution of t/t*, the fraction of a source's lifetime at which 
it is observed, is just a constant for a given sample, then D ~ D z /2. More 
generally, if the t/t* distribution of a sample does not evolve with redshift, then 
Doc D z . 

3. Samples 

In order to use powerful extended radio sources for cosmology, the compari- 
son sample must be as homogenous as possible. Thus, only radio galaxies with 
178MHz power P 178 (qo = 0) > 3h~ 2 X 10 26 W Hz" 1 sr" 1 are considered [11]. 
The power cut insures that backflows and bridge distortions are not significant, 
and that the lobes are supersonic [6,10,14]. Although this method might be ex- 
tended to radio-loud quasars, current observations favor the use of radio galaxies 



The radio galaxies used for the comparison sample are the 3CR radio galaxies 
with FRII morphology listed in Table 1 of McCarthy et al. [15], subject to 
the power constraint mentioned above. To calculate the radio power of these 
sources, fluxes and spectral indices from Spinrad et al. [16] are used. Cygnus 
A, which was not listed by McCarthy et al. due to its low galactic latitude, is 
included in this study since it is the closest high-power FRII radio galaxy, and 
is well-studied. 

To compute D* } values for vl, a^, and B are needed, and those extracted 
by Wellman & Daly [17,18] are used here. Wellman & Daly reanalyzed the 




(2.8) 




(2.9) 



[6]- 
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z range 



0.0-0.3 
0.3-0.6 
0.6-0.9 
0.9-1.2 
1.2-1.6 
1.6-2.0 



# sources 



3 

15 
27 
19 
10 

8 



(a) 

(kpc/h) 
33±5 
89±15 
63±7 
46±8 
41±11 
25±6 



(b) 

(kpc/h) 
66±14 
177±42 
126±18 
91±21 
83±30 
50±17 




98±17 
74±8 
57±10 
55±14 
36±8 



(d) 
(kpc/h) 
68±14 
197±47 
149±21 
114±27 
110±40 
71±24 




114±19 
90±9 
70±12 
69±17 
44±10 



(f) 
(kpc/h) 
72±13 
227±54 
179±25 
141±33 
137±50 
88 ±30 



TABLE 1. Average physical sizes by redshift bin for the sample described in §3: (a) 
Core-Lobe and (b) Lobe-Lobe for (O,A o ) = (1.0,0); (c) Core-Lobe and (d) Lobe-Lobe for 
(ft,A o ) = (0.1,0); (e) Core-Lobe and (f) Lobe-Lobe for (ft,A o ) = (0.1,0.9). 



published maps of 14 radio galaxies, from [8,9], to obtain consistent values for 
the parameters needed in eq. (2.8). These sources are only a subset of the 
comparison sample, and in the future, D*s for more of this sample should be 
obtained [11]. 



4. Calculations &: Results 

Redshift bins were chosen to be of similar size, and the average core-lobe 
and lobe-lobe sizes for the whole comparison sample are listed in Table 1 for 
(O o ,A o )=(1.0,0), (0.1,0), & (0.1,0.9). The average lobe-lobe sizes are about twice 
the average core-lobe sizes, as expected. Figure 1 shows the average core-lobe 
size for the sample. Except for the lowest redshift bin, the average size of these 
radio galaxies decreases with increasing redshift. 

Each redshift bin has an average physical size D. The data needed to calculate 
D* are only available for small number of sources, and the determination of 
typical D* for a bin based on one or two sources could be misleading [11]. In 
order to compare the evolution of D* to D } the ratio D / ' D* is calculated for each 
lobe and then each source, using D from the appropiate redshift bin. 

In principle, for the correct values of f3 and the cosmological parameters, the 
ratio D/D* is constant and independent of redshift. For example, the ratio 
D/D* assuming f3 = 2.0 and (O o ,A o )=(0.1,0) is shown in Figure 2 for each lobe 
with a calculated D*. The ratio in Figure 2 is nearly constant for all redshift, 
and similar results are found for different cosmological parameters [11]. 

At a given redshift, D/D* oc (a r) 3 ^ 7+2 ^ 3 [6], and for f3 > 1, the ratio should 
depend on (a r) to a power greater than one. Since the normalization of the ratio 
is arbitrary, and (a r) = H ~ /(0 , A,z), H does not enter in the analysis. 

Using the three sets of cosmological parameters mentioned above, and letting 
f3 vary, D/D* is fitted to a constant using the weighted mean, and a reduced 
X 2 is calculated. The values for the reduced \ 2 as a function of f3 including 
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FIGURE 1. The average core-lobe size, D, for the sample in Table 1. 
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Figure 2. D/D* for each lobe with (3 = 2.0 and (O o ,A o ) = (0.1,0). 



the Cygnus A point are shown in Figures 3 & 4 for the core-lobe and lobe-lobe 
calculations respectively. 

The best-fit value of f3 is about two in each of the three cosmologies studied. 
This is consistent with the result obtained by Daly [6,7], f3 = 1.5 ± 0.5, based 
on a different argument. In calculating the reduced \ 2 f° r the core-lobe com- 
parison, the two lobes of a given source are treated as statistically independent 
which is not, strictly speaking, correct, whereas the lobe-lobe comparison is valid 
statistically [11]. 

A low universe is favored over a flat, matter-dominated universe as can 
be seen in both Figures 3 & 4. In principle, this method is sensitive to the 
difference between a A-dominated or curvature-dominated cosmology, but both 
of these models fit the current data well [11]. 

Figures 5 & 6 are to be compared to Figures 3 & 4, and show the reduced 
X 2 's excluding Cygnus A. This eliminates the lowest redshift bin, which has the 
fewest sources to calculate D . These figures show that the results are insensitive 
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FIGURE 3. Reduced Chi-squared for Core-Lobe comparison including Cygnus A. 
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FIGURE 4. Reduced Chi-squared for Lobe-Lobe comparison including Cygnus A. 

to the lowest redshift bin. D for the lowest redshift bin and D* for Cygnus A are 
both low compared to the adjacent redshift bin, yet the ratios are comparable 
to those at higher redshift for f3 ~ 2 [11], as is evident in Figure 2. 
More details of this analysis and additional figures are given by [11]. 

5. Conclusions 

The results presented here indicate the promise of this method. The fact 
that Cygnus A seems to fit in the sample in terms of this analysis and other 
properties, such as Lj and n a [6,7] is encouraging. This analysis favors a low- 
density universe that is curvature-dominated or A-dominated over a flat, matter- 
dominated universe. 

Computing D* } s for more sources and increasing the size of the comparison 
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FIGURE 6. Reduced Chi-squared for Lobe-Lobe comparison excluding Cygnus A. 

sample are essential for the use of this method to estimate cosmological param- 
eters and to test the model for internal consistency. Thus, a long term goal is to 
expand the number and redshift range of both samples. 
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cussions. A special thanks goes to Greg Wellman for access to his data and 
for valuable discussions. This work was supported by an NSF Graduate Fellow- 
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